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Limitations  of  MD  simulation 


•  Discrepancy  between  MD  simulation  and  experimental 
results 

1.  Length  Scale 

-  Recent  billion  atoms  simulation  (lxlxl  jam  for  metal) 

-  Most  of  MD  simulation  is  nanometer  scale  (In  many  case,  using  only  a 
small  part  of  actual  specimen  for  simulations) 

2.  Time  scale 

-  Most  MD  simulation’ s  time  step  is  femto  second 

-  Current  time  step  extending  methods  can  increase  time  step  about 
factor  of  15 

•Necessity  to  increase  length  and  time  scale  of  MD  simulations 
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Method  and  framework 


•  Equivalent  crystal  lattice  method 

•  New  inter-atomic  potential  parameters 

-  Tersoff  potential  (T3) 

•  Comparison  of  properties 

-  Cl  1,  C12,  C44,  Bulk  modulus,  Total  lattice  energy. 

•  Time  step  calculations 

-  Maxwell-Boltzman  distrubution 

•  Simulation  details 

-  Construction  of  suprtcells,  NST  ensemble, 
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Equivalent  crystal  lattice  method 


4 


F=LwM“r-fr'rf 


i—l 


-  Lattice  parameters: 
5.43, 10.86, 108.6,  and 
1086  A 

-Number  of  atoms  in  the 
unit  cell  are  the  same 

—Fraction  coordinates  in 
the  unit  cell  are  the  same 

-Potential  parameters  are 
found  by  minimizing  F 


Nobs  is  the  number  of  observables,/"'  and/"7'  are  the  target  and  calculated  values  of 
the  observables,  and  wt  is  the  weight  factor  for  the  given  observables 
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Parameters  for  Tersoff  potentials 


Lattice 

Constant 

5. 43 A  (T3) 

10.86A 

108.6A 

1086A 

A  (eV) 

1830.8 

9911.8 

9682272.8 

12256245358.6 

B(eV) 

471.18 

2500 

2500000.0 

2846954605.4 

A  (A1) 

2.4799 

1.2125 

0.1195 

0.01235 

/I,  (A1) 

1.7322 

0.8166 

0.08071 

0.008282 

a 

0 

0 

0 

0 

n 

0.78734 

0.78734 

0.78734 

0.78734 

p 

1.0999E-06 

1.0999E-06 

1.0999E-06 

1.0999E-06 

n 

0.78734 

0.78734 

0.78734 

0.78734 

c 

100390 

100390 

100390 

100390 

d 

16.218 

16.218 

16.218 

16.218 

h 

-0.59826 

-0.59826 

-0.59826 

-0.59826 

'  My 

Ly  =/c <  rij  )[a,y/jl  ( 'Vy  )  +  ^ijf  A  <  'i]  >  1 
fR{r)=  A  exp(  —  Xjr )  , 

/^(r)=  —  5  exp(  —  A2r)  , 


/c(r): 


1,  r<R-D 

1 T 


T~Tsin 


(/■-*)//> 


,  i?  —  Z>  <r  <R  +D 


0,  r>jR+Z) 


(lc) 


fc,7  =  (l+|8"fj)-1/2"  , 

£</  =  2  /c(r,k)*(®yt)exp[A.5(ry-rft)3]  ,  (Id) 

*  (gti,y) 

g(d)=l+c2/d2  —  c2/[d2  +  (h  —  cos0)2]  , 

aij={\+anrj'-j)~W2n  , 

,  (le) 
Vij  =  2  /c  ( rik  )  exP[  -  ■ rik  )■ 1  ]  . 

*  M*.y> 
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Results  comparison 


( 


Properties 

Target  values 
(experimental) 

Tersoff  (T3) 

New  potential 

Lattice  parameters(A) 

5.43  A 

10.86  A 

108.6  A 

1086  A 

cu  (GPa) 

165.7 

142.5  (14) 

137.0  (17) 

134.86  (18.6) 

136.6(17.5) 

C12  (GPa) 

63.9 

75.4  (18) 

70.0  (9.7) 

68.13  (6.63) 

73.5  (15) 

C44  (GPa) 

79.6 

69  (13.3) 

67.42  (15) 

67  (15.8) 

64.8  (18.6) 

B  (GPa) 

98 

98  (0.0) 

98.5  (0.6) 

90.37(7.7) 

94.5  (3.5) 

Lattice  E  (ev) 

-37.36  (5.43A) 
-298.88  (10.86A) 
-298880  (108.6A) 
-298880000 

(1086A) 

-37.04(0.9) 

-289.4  (0.14) 

298387  (0.1) 

-298935530(0.05) 

Atomic  mass 

28.0855 

28.0855 

224.684 

224684 

224684000 
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Silicon  nanowires  Simulation  detail 


Supercell 

size 

5.43  A 

10.86  A 

108.6  A 

1086  A 

4x4x32 

2.17x2.17x17.36  nm 
(4096) 

4.34x4.34x34.8  nm 
(4096) 

43x43x348  nm 
(4096) 

0.43x0.43x3.48pm 

(4096) 

8x8x32 

4.34x4.34x17.36  nm 
(16384) 

8.7x8.7x34.8  nm 
(16384) 

87x87x348  nm 
(16384) 

0.87x0.87x3.48pm 

(16384) 

16x16x32 

8.7x8.7x17.36  nm 
(65536) 

17.4x17.4x34.8  nm 
(65536) 

174x174x348  nm 
(65536) 

1.74xl.74x3.48pm 

(65536) 

-  Applying  tensile  load  constant  strain  rates: 

0.01  %/step  and  0.001  %/step 

-  NST  ensemble  (Constant  particles,  temperature,  stress) 


-  Temperature  :  300K,  Nose-Hoover  thermostat 
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Simulation  results  between  5.43  and  10.86A  equivalent  lattice 


40 


CL 

o 

C/3 

CO 

CD 

CO 


5.43A  8x8x32  supercell 

40 

co  30 
CL 

CD 

5.43A  8x8x32  supercell 
10.86A  4x4x16  supercell  / 

^20 

CO 

L^v/ 

1 

^^M0.86A  4x4x16  supercelj 

\  in  9n  a 

-S  io 

80  . 
n  0 

Strain  (%) 


0 


Strain  (%) 


40 


(a)  strain  rate  0.01  %/step  (b)  strain  rate  0.001  %/step 

Dimension:  4.34  x  4.34  x  17.36  nm 
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Fracture  results  between  5.43  and  10.86A  equivalent  lattice 


8x8x32  at  32% 


4x4x16  at  30% 


8x8x32  at  34% 


4x4x16  at  32% 


(a) 


5.43  A 


10.86 


(b) 


Dimension:  4.34  x  4.34  x  17.36  nm 
strain  rate  0.001  %/step 
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Statistical  analysis:  Coefficient  of  correlation  (R) 


SSR 

SSTO 


rw 

V  SSTO 


ZA  —  z 

(Yi  -  Y )  :  Regression  sum  of  square, 

-  A 

where  Y :  mean  value,  Yt :  value  from  the  fitted  line 

SSE  =z_l(Yi-Yi)  :  Error  sum  of  square 
SSTO  =  SSE  +  SSR  :  Total  sum  of  square 


Properties 

Predictor  variables 

Diameter 

Surface 

area 

Volume 

Surface  to 
volume  ratio 

Young’s 

modulus 

0.14 

0.10 

0.15 

-0.74 

Fracture 

stress 

0.30 

0.25 

0.26 

-0.85 

-  The  surface  to  volume  ratio  is  most  critical  factor 
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Young's  Modulus,  [GPa] 


Young’s  modulus  transition  diameter  (Si  nanowire) 


100 


70 


'  r'X'  "I 

■ 

E  !  ! 

proposed 

™  i  Current 

Approach 

t  !  MD 

and 

■  ■  Simulation  ■ 

experiments 

r  1  1 

■ 

.Li  L 

- ✓/ - 

0  30 


100 

Diameter  (nm) 


2000 


Discrepancy  b/w  current  MD 
and  experiments 

-Current  MD  predicts  4-30nm  for 
transition  diameter  (Park  et  al., 
MRS  bulletin,  34,  2009) 

-Experiments  measured  lOOnm 
for  transition  diameter  (Li  et  al. 
Appl.Phys.Lett,  83,2003  qnd  Han 
et  al.  adv.  Mater.,  19,  2007) 

-Our  simulation  predicted  around 
lOOnm  for  transition  diameter 
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Polysilicon  simulation  detail  (Number  of  atoms  and  sizes) 


5.43  A 

10.86  A 

108.6  A 

1086  A 

8  grains 

49356 

(101x101x101  A) 

49356 

(202x202x202  A) 

49356 

(2020x2020x2020  A) 

49356 

(20200x20200x20200  A) 

64  grains 

48058 

(101x101x101  A) 

48058 

(202x202x202  A) 

48058 

(2020x2020x2020  A) 

48058 

(20200x20200x20200  A) 

-  Constructed  bulk  polysilicon  and  surfaced  polysilicon 


-  Applying  tensile  load  with  constant  strain  rates: 
0.01  %/step  and  0.001  %/step 


-  NST  ensemble  (Constant  particles,  temperature,  stress) 


-  Temperature  :  300K,  Nose-Hoover  thermostat 
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Number  of  grain  and  surface  effect  on  properties 


(a) 


IE-3  0_01 

Surface  to  Volume  ratio 


8. 
o 

“25H 
(0 
CO 
0) 

^  20-| 
0) 

3 

o  15J 
co 


(b) 


8  grains 

7':x\ 


10- 
1E-4 


64  grains 


TT|- 

0.01 


IE-3 

Surface  to  Volume  rafio 


0.1 


(C) 


(d) 


(a) Young’s  modulus  with  strain  rate 

0.01  %/step 

(b)  Young’s  modulus  with  strain  rate 

0.001  %/step 

(c)  fracture  stress  with  strain  rate 

0.01  %/step 

(d)  fracture  stress  with  strain  rate 

0.001  %/step 


-  With  the  same  surface  to  volume  ratio  (same  dimension),  number  of 
grain  changes  the  Young’s  modulus:  Grain  size  effect 
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Comparison  with  experiments 


-Koskinen  et  al.  measured 
fracture  stress  with  various 
grain  size 

-  Inverse  Hall-Petch  mechanism 


-  Our  results  agree  with 
Koskinen’ s  experiments 

0  20  40  60  80  400  500  600  700 
Grain  size  (nm) 


CO 

CL 

CD 

in 

w 

<D 

w 

o> 


o 

co 


8 

6 

4 

2 

n. 


m=0.0474 

■ - 

Experiment 
by  Koskinen 

L- 

Most  ■ - 

experiments 


m=0.079 


/ 


m=0.112 

U _ 


-  Our  fracture  stresses  are  converted  assuming  typical  experimental  strain  rate  2E-3/sec 


a 


(  •  \m 

£ 


a 


0 


V^o  J 


Q  :  stress,  £  :  strain  rate,  m  :  strain  rate  sensitivity  index 
(  0.0474  <  m  <  0.1 12,  avgerge  m  =  0.079) 
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Limit  of  surface  and  grain  size  effects  on  polysilicon 


strain  rate  0.01  %/step 


-  Grain  size  effect  and  surface  effect  is  decreasing  and  converge  to  0  !! 


-  About  1  |Lim  of  the  grain  size  and  23  pm3  of  specimen  size  are  the  limit 
above  which  the  Young’s  modulus  and  fracture  stress  is  not  influenced 
by  grain  size  as  well  as  surface  for  polysilicon 

-  Agree  with  experiments  by  Fancher  et  al. ,  J.  Mater.  Sci.,  36,  2001, 
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Fracture  transition  from  brittle  to  ductile  with  surfaces 


(a) 
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(a)  bulk  polysilicon 

(b)  surfaced  polysilicon 


(b) 


101  x  101  x  101  A  8  grain  polysilicon 
with  strain  rate  0.001  %/step 

-Surfaces  make  materials 
more  ductile  when  compared 
to  bulk  materials 

-Slip  induced  fracture  is 
observed  inside  grain 
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Fracture  behaviors  for  equivalent  lattices 


( 


(a) 


(a)  5.43  A  lattice 

(b)  10.86  A  equivalent  lattice 

o 

(c)  108.6  A  equivalent  lattice 

(d)  1086  A  equivalent  lattice 

64  grain  bulk  polysilicon  with 
strain  rate  0.001  %/step 


-Fracture  initiations  are 
identical  for  all  lattices 


-Fracture  initiation, 
propagation,  final  pattern  are 
identical  b/w  5.43  Aand  10.86  A 


<C) 
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-Temperature  s  are  higher  than 
300K  for  108.6A  and  1086A 
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Intergranula  vs.  transgranula  depending  on  strain  rate 


Sz2  gslsrl  and  sr2  - 
surface 

202  x  202  x  202  A  8  grain 
surfaced  polysilicon  with 

o 

strain  rate  with  10.86A 
equivalent  lattice 
(a)  0.01  %/step  (b) 

0.001  %/step 


UNIVERSITY  OF 

NOTRE  DAME 


Aerospace 

Mechanical 

Engineering 


(±1 

Multiphysics  Labs 


-Fracture  propagation 
depends  on  strain  rate 

-Mostly,  intergranula 

-Transgranula 

observed 

Notre  Dame  (http://wwwjid.edu/~vtomar) 


Conclusion 


•  By  using  equivalent  crystal  lattice  method, 

•  Time  step  increased  about  factor  of  500 
(From  0.001  to  0.5fs) 

•  Length  scale  increased  factor  of  200 
(5.43  A  to  1086  A) 

•  The  surface  to  volume  ratio  is  the  most  important  factor  at 
nanometer  and  micron  scale 

•  MD  simulation  results  are  more  close  to  experiments  if 
actual  specimen  dimension  is  used  for  MD  simulations 

-  Young’s  modulus  transition  diameter  (Si  nanowire) 
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Conclusion 


•  MD  simulations  of  polysilicon  with  grain  sizes  ranging  from 
3.4nm  to  1.3  pm  are  made  possible  using  equivalent  lattices 

•Inverse  Hall-Petch  mechanism  is  observed  for  polysilicon 

-  Agrees  with  experiments,  opposite  to  metals  and  alloys 

•  The  Young’s  modulus  is  not  influenced  by  grain  size  when 
grain  size  is  larger  than  7nm  for  bulk  poly  silicon 

-  Verified  by  DFT  calculations 

•  The  Young  modulus  is  strongly  influenced  by  grain  size  for 
surfaced  polysilicon 


•The  existence  of  surfaces  make  materials  more  ductile  when 
compared  to  bulk  materials 
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